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ABSTRACT 

High- velocity clouds (HVC) , fast-moving ionized and neutral gas clouds found at high 
galactic latitudes, may play an important role in the evolution of the Milky Way. The 
extent of this role depends sensitively on their distances and total sky covering factor. 
We search for HVC absorption in Hubble Space Telescope high resolution ultraviolet 
spectra of a carefully selected sample of 133 AGN using a range of atomic species in 
different ionization stages (e.g., O I, C II, C IV, Si II, Si III, Si IV). This allows us 
to identify neutral, weakly ionized, or highly ionized HVCs over several decades in 
H I column densities. The sky covering factor of UV-selected HVCs with Ii^lsrI ^ 90 
kms~^ is 68% ± 4% for the entire Galactic sky. About 74% of the HVC directions 



have iV(H l) < 3 x lO^^ cm'^ and 46% have 7V(H l) < 



10^ 



We show 



that our survey is essentially complete, i.e., an undetected population of HVCs with 
extremely low H (H I+H II) column density is unlikely to be important for the HVC 
mass budget. We confirm that the predominantly ionized HVCs contain at least as 
much mass as the traditional H I HVCs and show that large H I HVC complexes have 
generally ionized envelopes extending far from the H I contours. There are also large 
regions of the Galactic sky that are covered with ionized high-velocity gas with little 
H I emission nearby. We show that the covering factors of HVCs with 90 ^ \vi^sr\ ^ 

170 kms~ drawn from the AGN and stellar samples are similar. This confirms that 
these HVCs are within 5-15 kpc of the sun. The covering factor of these HVCs drops 
with decreasing vertical height, which is consistent with HVCs being decelerated or 
disrupted as they fall to the Milky Way disk. The HVCs with |wlsr| ^ 170 kms^^ 
are largely associated with the Magellanic Stream at b < 0° and its leading arm at 
6 > 0° as well as other large known H I complexes. Therefore there is no evidence in 
the Local Group that any galaxy shows a population of HVCs extending much farther 
away than 50 kpc from its host, except possibly for those tracing remnants of galaxy 
interaction. 

Key words: Galaxy: halo - Galaxy: evolution - Galaxy: structure 



1 INTRODUCTION 

The ability of galaxies to form stars depends sensitively on 
the content and physical conditions of their gas, which in 
turn is dictated by internal effects, such as feedback, as 
well as by the external interaction of galaxies with their 



* Based on observations made with the NASA/ESA Hubble 
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tute, which is operated by the Association of Universities for Re- 
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surroundings. A galaxy may exchange mass with its circum- 
galactic environment through the infall of intergalactic mat- 
ter or of remnants resulting from galaxy interactions, and 
through the outflows driven by stellar and AGN feedback. 
Despite many theoretical and observational efi'orts (a non 
exhaustive list of references includes |Oort||1970| [Steidel fc 



2007 



Peek 



lSargent"1992- 'Keres ot al."2005" "Bou che et al 

et al. 2008; Dckcl ct al. 2009; Prochas kafc Hennawi| [2009| 



Oppenheimer et al.||2010{ 



jet al.||2011[ [Tripp et aL 



Fumagalli ot al. 2011, Marinacci| 
; |Tumlinson et al. 1,2011,) , the 



2011 



jogf 

acci| 



cycle of baryons in and out of galaxies and its implications 
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for galaxy evolution are not yet fully understood. This ex- 
change of baryons between galaxies and their environments 
must, however, result in a net gain of mass if galaxies are to 
form stars over many billions of years. 

In our own Milky Way, it has been long known that 
there must be streams or flows of gas through the Galac- 
tic halo. Signatures of these flows are revealed by the high- 
velocity clouds (HVCs), clouds moving in th e local standard 
of rest (LSR) frame at |wlsr| ^ 90 kms~^ (Wakker & van 
Woerderi||1997 l. The unknown distances to these clouds for 
a long time meant we did not know if the HVCs were associ- 
ated with streams occurring near the Milky Way (e.g. 
1970[) rather than the IGM of the Local Group (e.g. 



Oort 



Blitz 



et al.||1999||Gnat fc Sternberg||2004l ). Furthermore, the dis- 
tances to the HVCs are required for quantifying their basic 
physical properties (several scale with the distance, e.g., the 
mass M oc dP). Since the last large summary on the distances 
and metallicities of the H I HVCs published 10 years ago 
( Wakker||2001 1 , our knowledge of these clouds has dramat- 
ically improved. Distances of several individual large HVC 
complexes seen in H I emission are now determined, plac- 
ing them at heliocentric distances of about 4-15 kpc (e.g.. 



Ryans et aLl |T997; Wakkor 2001) |Thom e t aL 2006 
Wakker et al 



2008 



2007j _2008_, Lehner & Howk:2010_; Smoker et_ 



al.|2011 1 . Among the large complexes (in terms of mass and 
solid angle), only the Magellanic Stream (e.g., Putman et al. 
19981 [Briins et al.|2005| [Nidever et al.|2008| ) is much farther 
away, possibly extending to 80-200 kpc according to recent 
numerical models of the interaction of the Large and Small 
Magellanic Clouds (IBesla et al.||2012l). 



With the advent of sensitive high resolution UV spectro- 
graphs including the Far Ultraviolet Spectroscopic Explorer 
(FUSE) and those of the Hubble Space Telescope (HST), it 
has also become obvious that HVCs are both neutral and 
ionized in view of the high-velocity absorption detected in 
resonance lines of neutral, weakly ionized, and highly ion- 
ized species in the spectra of AGN (e.g.. 



et al.||2003| .Wakker et al.| 2003 Collins et al. 2005 |Fox et 



Sembach et al. 



1995 



1999 , 2000; Sava ge et al.|2000| [Lehner et al.|2001| [Sembach 



al.|20b6[[ShuU et al.|2009[(Richter et al.|2009[ [Tripp fc Song 



20121. The importance of ionized gas can be demonstrated 



simply by considering the fraction of the sky covered by 
HVCs at H I column sensitivities. While only 18% is cov- 
ered at a sensitivity of logiV(H I) > 18.5 (|Wakker||1991|), 
a limit easily detected in H I 21 cm observations, deeper 
H I emission observations (log A*'(H I) J> 17.9) yield a factor 
larger covering factor, 37% (Murphy et al._1995^,Lockman et 
aL]|2002 1 . A much larger fraction of the sky is even covered 
by the HVCs seen in absorption, which are sensitive to much 
lower Af(H I): for the HVCs with Q VI absorption, the HVC 
covering factor is about 60-75% ( Sembach et al.||2003 Fox 
et al. 20061, and a high covering factor was also found for 



the HVCs detected via the strong line of Si III ( Shull et al, 
2009{|Comns et al.|2009] hereafter CSG09) 



Many of these HVCs observed in absorption have been 
often dubbed as ionized HVCs or highly ionized HVCs since 
in many cases a large fraction of the gas is photoionized or 
coUisionally ionized, i.e., A'^(H II) ^ A''(H I). In this work, 
we simply label them as HVCs. To differentiate them from 
the HVCs observed in H I emission, we refer to the latter 
category as H I HVCs. We emphasize that the only real 
distinction with this definition is that the H I HVCs are 



detectable via 21-cm emission. Thus they have a clear H I 
sensitivity cut-off determined by the depth of the observa- 
tions and limitation of radio emission studies. On the other 
hand, for the HVCs (i.e., those seen in UV absorption), the 
amount of H I can vary greatly as demonstrated by the ob- 
served range in the H I column density; A''(H I) <: 10^^ 
cm~^ to A'^(H I) > 10^° cm~^ derived from Lyman series 



absorption (Fox et al. 2006 Zech et al. 20081 and H I 21 



cm emission observations toward the AGN; virtually any 
N{¥^ I) column densities are observed in the UV-selected 
HVC. However, even if the H I column density can be small, 
the total H (H I-fH II) column density of the HVCs can 
stiU be substantial, > lO^'^-lO^^ {Z/Zq)~'^ c m~^ on average 
(e.g., [Sembach et aL)2003| [Shull et al.|2009[ ). As for the HI 
HVCs, direct distance constraints are required for determin- 
ing the masses of the HVCs seen in UV absorption and for 
characterizing their role in the evolution of the Milky Way 
and their relationship with the H I HVC complexes. If they 
are in the Milky Way within several kpc, the HVCs seen in 
UV absorption are not only more important in number and 
mass, but their relative proximity may allow them to de- 
scend into the disk to fuel new generations of star formation 
(e.g., [ShuU et al.|2009| [Lehner fc Howk[[2011 l. 



[Lehner & Howk[(|2011| hereafter LHll) ruled out whole 



classes of HVCs using a novel and powerful method to statis- 
tically constrain the distance of the entire HVC population 
rather than just a few individual clouds. This was done by 
comparing the UV absorption line detection rates of HVCs 
in COS and STIS spectra of stars (distance limited) and 
AGN (no distance limit). They selected 28 stars (hereafter, 
the LHll stellar sample) to be at the largest practical dis- 
tances in order to distinguish between Milky Way halo and 
Local Group origins for the HVCs, resulting in targets at 
\z\ > 3 kpc from the Galactic plane and heliocentric dis- 
tances 3 ^ d <; 32 kpc (although most of the stars are within 
(d) — 11.6 ± 6.9 kpc). The stellar sightlines were carefully 
selected to be well distributed across the high latitude Galac- 
tic sky and avoid biases that a priori favor or disfavor the 
presence of HVCs along each sightline. They identified HVC 
absorption in the stellar spectra using the strongest lines of 
atomic, singly, and doubly ionized species available in the 
UV bandwidth (e.g., C II A1334, C IV A1548, Si II A1260, 
Si III A1206, O I A1302), most showing high [Si II/O I] ratio, 
consistent with these HVCs being mostly ionized. The HVC 
sky covering factors of the HVCs with 90 ^ I^ilsrI < 170 
kms~^ determined from the stellar and AGN samples are 
fc — 0.50 and <. 0.65, respectively (LHll). Unfortunately 
the HVC covering covering factor for the AGN sample suf- 
fers from poorly understood selection criteria, resulting only 
in an upper limit on fc- Nevertheless the similarity in the 
covering factors between the distance limited sample and 
the sample with no distance limit imply that most of the 
HVCs with 90 ^ JulsrI ^ 170 kms~^ are within heliocen- 
tric distances 5-15 kpc (no HVC at |wlsr| > 200 kms~^ was 
found in the stellar sample). Given the importance of HVC 
distances and the importance of the covering factor in deter- 
mining d, one of our goals here is to robustly determine the 
covering factor of the HVCs from a well understood sample 
of AGN that can be directly compared to the LHll stellar 
sample. To search for the HVCs, we do not focus on a single 
ion, but on a suite of atoms and ions (e.g., O I, C II, C IV, 
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Si II, Si III, Si IV, Al II) following the same method used 
by LHll. This allows us not only to securely detect HVC 
absorption, but also, as we will show, to sensitively search 
for HVCs at all levels of ionization. 

Here we present the results of our new study of HVCs, 
with the following organization. In ^we describe our sam- 
ples of AGN, the data and their sensitivity, and the method 
for determining the presence of HVCs. In 33] we present our 
findings for the covering factors of the HVCs and VHVCs. 
With this new sample of AGN, we strengthen the LHll re- 
sults, and we extend their results by using stars at smaller 
distances to study the dependence between the covering fac- 
tor and the distance. With our much larger sample of AGN 
and much better sky coverage than the sample used in LHll, 
we are also able to gain insight into the origins and distances 
of the very high-velocity clouds (VHVCs, |i;lsr| > 170 
kms~^). In §4|we connect the HVCs detected in H I emission 
and UV absorption based on their closed proximity on the 
sky, in velocity, and in space (distance). We discuss some of 
the implications of our results for the mass of the HVCs and 
their origins and fate in !|5] Finally, we present a summary 
in^ 



2 AGN SAMPLE AND DATA 

2.1 Description of the main AGN sample 

The 67 AGN in our main sample were drawn from 3 
HST/COS large programs: 11598 ("How Galaxies Acquire 
their Gas: A Map of Multiphase Accretion and Feedback in 
Gaseous Galaxy Halos", PI: Tumlinson), 11741 ("Probing 
Warm-Hot Intergalactic Gas at 0.5 < z < 1.3 with a Blind 
Survey for O VI, Ne VIII, Mg X, and Si XII Absorption 
Systems, PI: Tripp), and 12248 ("How Dwarf Galaxies Got 
That Way: Mapping Multiphase Gaseous Halos and Galac- 
tic Winds Below L*", PI: Tumlinson) Pj The aims of these 3 
programs are to target science beyond the Milky Way halo. 
Therefore the choice of the targets was completely random 
with respect to the locations of the H I HVCs. Hereafter we 
refer to this sample as the TT (Tripp-Tumlinson) sample. 

Information on the design and performance of COS can 
be found in Green et al. (20121. The data reduction and co- 



addition procedures are described in Thom et al. (2011 1 and 



Meiring et al. ( 2011 1, and we refer the reader to these papers 
for more details on data processing. Notably the individ- 
ual exposures were co-added in photon counts, rather than 
flux calibrated exposures. The shifts between the individual 
exposures were determined using the Galactic interstellar 
absorption lines. The exposures were shifted to a common 
reference and co-added. All the data were then shifted to 
the local standard of rest (LSR) frame. A visual comparison 
between the Galactic H I emission spectra from the Lei- 
den/Argentine/Bonn (LAB) survey ( Kalberla et al. 2005 1 



and the Galactic absorption indicates that the spectra were 
correctly shifted in the absolute LSR frame. The COS res- 
olution of the G130M and G160M gratings is R ^ 17, 000, 



which is adequate for identifying HVCs (Lehner & Howk 
2010 LHll). The typical signal-to-noise of the co-added 



spectra are about 4-9, > 20, and 9-13 per resolution ele- 
ment for programs 11598, 11741, and 12248, respectively. A 
few sightlines were not considered in the programs 11598 
and 12248 owing to the poor signal-to-noise of their spec- 
tra {<. 3) or excessive contamination from higher redshift 
absorbers. 

In Fig. Ill we show the sky distribution of the 67 TT tar- 
gets, and in Table[l]we list their Galactic coordinates. Most 
of the targets are situated in the northern Galactic hemi- 
sphere as the two largest samples (11598 and 12248) were 
selected from the Sloan Digital Sky Survey (SDSS). There 
is an excellent coverage of the Galactic sky in longitude and 
b > 30° . There is also good overall agreement in the sky cov- 
erage between the TT sample and sample of stars assembled 
by LHll as demonstrated in Fig. [l] 



2.2 Additional sample of AGN 

In order to achieve the best sky coverage possible, we also 
consider two additional samples: 1) 18 publicly available 
COS AGN spectra from the Guaranteed Time Observer 
(GTO) IGM program, and 2) 49 AGNs drawn from the STIS 
E140M and G140M AGN sample from CSG09. 

The GTO COS spectra were reduced following the pro- 
cedure described in §2.1| We only considered AGN that 
were initially targeted for IGM science, not the AGN from 
the GTO Galactic-HVC program to avoid obvious selection 
bias. The GTO sample is summarized in Table [2] 

The only explicit two criteria for selecting sightlines in 
CSG09 were that 1) S/N ^ 3, and 2) little apparent con- 
tamination from the QSO absorbers. However, we need to 
remove any potential bias from their sample, as, e.g., some 
of their AGN were initially targeted for HVC science (with 
a priori knowledge from past observations that an HVC was 
along the line of sight). In order to identify AGN originally 
targeted for HVC science, we reviewed the MAST program 
abstract and title for each target in CSG09 and found that 
7 sightlines were explicitly selected to study HVCs. Those 
were not considered further. From the appendix in CSG09, 
we determined that 6 lines of sight had an HVC detection 
solely based on Si III absorption (for the other HVCs, they 
have confirmation of the absorption using several species 
and/or transitions, i.e., consistent with our criteria from 
f 2.3 1. This differs from our method in which we only classify 



an HVC detection based on its detected absorption in more 
than one species. Considering these sightlines but classifying 
them with no HVC detection or removing these sightlines 
from the counting would not change the results apprecia- 
bly. We include these sightlines but do not count them as 



detections (see next in f2.3l 



^ The program 12248 is not completed yet; only data taken prior 
to November 2011 arc in this paper. 



2.3 Search for HVC absorption 

In order to search for the high-velocity interstellar absorp- 
tion (90 s; |i;lsr| s; 400 kms"^) in the TT and GTO sam- 
ples, we first used the strongest lines of atomic, singly, and 
doubly ionized species available in the COS UV bandwidth: 
C II A1334 (/A = 171), Si II A1260 (/A = 1487), and Si III 
A1206 (/A = 1967). The smaller /A for C II is compensated 
by the fact that C is about 8 times more abundant than Si in 
a standard solar abundance pattern. These ions have been 
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Figure 1. Aitoflf projection map of the survey directions for the TT sample. Sightlines are plotted in Galactic coordinates with longitude 
increasing from right to left. A colored circle indicates an HVC along the line of sight while open circle implies no HVC along the sightline. 
Circles with a cross indicate that the sightline is part of the high-sensitivity sample with 10 ^ Wiim Sj 42 mA. The velocity value is 
color coded following the horizontal color bar. The star symbol shows the positions of the stars from the sample of LHll where a filled 
symbol indicates that an HVC is detected in the foreground of the stars, while an open symbol shows the absence of HVC along the 
stellar sightline. 



shown to be extremely powerful for finding low H I column 
density HVCs, and more precisely HVCs with virtually any 
H I column densities (see below). The use of multiple tran- 
sitions is crucial when using high redshift AGN because the 
likelihood of unrelated QSO absorber contamination is not 
negligible. Our rule is that an HVC absorption is defined as 
such only if it is seen in absorption in at least two species 
or atomic transitions. The same rule was used by LHll for 
their stellar sample. To confirm the detection or to help de- 
termine the velocity of the HVCs, we also used tracers of 
neutral gas (e.g., O I), weaker transitions of Si H, and trac- 
ers of more highly ionized gas (e.g., C IV). The strong line of 
the C IV doublet has /A = 147, comparable to the strength 
of C II, but C IV is more suitable than C II if the gas is 
highly ionized. In Fig [2] we show four sightlines from the 
TT sample, highlighting the diversity of UV-selected HVCs. 
The absorption profiles in this figure were normalized using 
Legendre polynomials with degrees d Sj 3. 

Fig. [2] shows that specific examples that the ionization 
conditions change fro sightline to sightline with in some cases 
only the detection of C II, siii, and Si HI (no high ion), and 
in other cases only detection of Si HI and C IV (no low ion). 
To illustrate this change in the ionization conditions, we 
show in Fig.^a Cloudy (Ferland et al. 1998 1 photoionization 
model that predicts the column densities of H I, C II, C IV, 
Si II, Si III, and Si IV as a function of ionization parameter 
{U = n-y/n-a, where n~f and nn are the densities of ioniz- 
ing photons and hydrogen). The radiation field is from the 
combined Milky Way and UV background radiation fields 



at a Galactocentric radius {Rg) of 10 kpc and a a-height of 
5 kpc ( Fox et aL]|2005 20101. These distances were chosen 
because they are close to the observed values of the HVCs 
(this paper, LHll 



Thom et al. 2006, 2008 



Wakker et al. 
[2008, 2007). We consider a total H column density of 10" 



cm and a metallicity of a 1/3 solar (see f 5.1 1. It is evident 



from this figure that the combination of several ions of the 
same and different species is extremely powerful for finding 
HVCs; e.g., as U increases, the amount of C II diminishes 
and may become undetectable, but then C IV can be de- 
tected. Therefore the available UV diagnostics allow us to 
search sensitively for any types of HVCs, neutral, weakly, or 
highly ionized, and this is well illustrated in Fig. [2] Also note 
that the H I column density decreases by a large amount 
from about lO^'' to lO^'^ cm'^ for -5 ^ log?/ ^ -1. The 
total H column density is by definition always 10^^ cm~'^ in 
this illustration. 



The fourth column of Table [T] summarizes the detec- 
tion or non-detection of HVCs along the sight lines of the 
TT sample, including the velocity of detected HVC ab- 
sorption. The velocities listed in this table were estimated 
from the apparent optical depth-weighted mean velocity 
(v) = j vTa{v)dv/ J Ta{v)dv. A colou iu this column high- 
lights that the velocity is severely blended with the lower 
velocity gas. A velocity with the superscript c indicates that 
the HVC was only detected in Si HI and C IV; there are only 
5 such cases. 
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Figure 2. Examples from the TT sample showing the variety of HVCs and S/N. The dotted vertical lines show the HVC detections. 
The only sightline that passes through an H I HVC complex (A) is J0950+4831, but only the negative component is associated with 
complex A along this direction. There is no H I HVC complex near the other directions (however, while complex L is several degrees 
away from J1437— 0147, it has similar LSR velocities as seen toward this sightline). Also note how the ionization conditions can change 
with, e.g., in some cases absorption only in singly and doubly ionized species (J0751+2919) and in other cases only in doubly and triply 
ionized species (J1233— 0031). 




near C H or Si H would be very similar (although the actual 
value of Wiim of C II or Si II may be different from those 
listed in Table [I] for a given sightline). To estimate Wiim, we 
measured the equivalent width (and la error) over the same 
velocity range (~ 80 kms^^) that is usually observed when 
an HVC is detected. The 3a upper limit on the equivalent 
width is defined as the la error times three. The mean and 
median of Wum are 42 mA, and the dispersion around of 
the mean is ad ~ 16 mA. Assuming the absorption line lies 
on the linear part of the curve of the growth, (W^um) cor- 
responds to logiV(C HA1334) ~ 13.3, logiV(C IVA1548) ~ 
13.0, logiV(SiHA1260) ~ 12.4, logiV(Si HI A1206) ~ 12.3. 



Figure 3. Cloudy simulations predicting the column density of 
the observed ions against the ionization parameter U = nry/nn. 
The radiation field is from the combined. Milky Way plus UV 
background radiation field at the position Rq = 10 kpc and 2 = 5 
kpc. Note that in this Cloudy simulation Af (H) is fixed and Af (H I) 



2.4 Sensitivity 

Knowing the sensitivity of our sample is key for determining 
if the detection of the HVC is dependent on the S/N of the 
spectra and for comparing or combining our sample with 
other samples. Below we demonstrate that the TT sample 
is sensitive enough to detect HVC absorption at a level pre- 
viously detected in other samples. 



2.4- i Sensitivity of the TT sample 

In Table [T] we report the 3a upper limits on the equivalent 
width (Wiim) obtainable for each line of sight from the TT 
sample and display its distribution in Fig. [4] We estimated 
Wiim near the Si III absorption, but the distribution of Wiim 



2.4-2 Comparison with previous samples of detected HVCs 

In order to compare our results with others, we estimated 
the equivalent widths W\ of C II for the HVC component 
when it is detected. We chose C II because this ion is cleanly 
observed in the stellar spectra. For example, in stellar spec- 
tra. Si III A1206 cannot always be estimated owing to con- 
tamination by the star itself or because the blended inter- 
stellar and stellar Lya absorption can severely depress the 
flux near Si HI. We did not attempt to separate severely 
blended HVC components and only measured W\ when the 
HVC absorption was not entirely blended with the lower 
velocity absorption. In Fig. |4] we show the distribution of 
Wx{C II) in the TT sample. The strength of most of the de- 
tected HVC absorption features is above (Wiim) + 2(7^ and 
is strictly larger than (Wiim). 

Using the stellar sample of LHll, we similarly estimated 
W\{C II) when a HVC is detected. Although the average 
S/N of the stellar spectra is higher than the S/N of the 
AGN spectra in the TT sample, the equivalent widths of 
HVCs in the stellar sample are such that Wx > (VKum) + 2(T(J 
(see Fig. E|. This implies that the detection rates of the 
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Table 1. The AGN sample 



Table 2. The GTO IGM sample 



Name 


I 


b 

(°) 


„HVC 
(kms-1) 


(mA) 


pt, 


J1342+1844 


0.13 


+75.48 




47.7 


2 


J1524+0958 


14.89 


+50.12 


-106 


19.1 


3 


J1409+2618 


34.67 


+ 72.59 




10.0 


3 


J1451+2709 


39.62 


+63.43 


-95 : 


41.4 


2 


J1330+2813 


42.37 


+81.23 




42.1 


1 


J1619+3342 


54.59 


+45.18 




30.7 


1 


J1445+3428 


56.74 


+64.59 


-101 


47.9 


1 


J1553+3548 


57.26 


+50.67 




40.4 


1 


J1555+3628 


58.32 


+50.27 




49.7 


1 


J1632+3737 


60.34 


+42.94 


-152 


9.5 


3 


J1330+3119 


61.27 


+80.41 




36.7 


2 


J1435+3604 


61.52 


+66.26 




60.2 


1 


J1550+4001 


63.96 


+50.94 


-108 


81.1 


1 


J1419+4207 


78.58 


+66.66 




52.8 


1 


J1342+3829 


83.11 


+ 74.40 


-102 


27.3 


2 


J2257+1340 


85.28 


-40.73 


-370 


64.6 


1 


J2345-0059 


88.79 


-59.39 


-276 


48.5 


1 


J1341+4123 


90.59 


+ 72.48 


-93 


27.0 


3 


J1322+4645 


107.71 


+69.44 


-100 : 


47.0 


1 


J0042-1037 


115.13 


-73.36 


-98 


53.4 


1 


J1241+5721 


125.48 


+59.73 


-158, -113 


37.0 


1 


J1233+4758 


131.24 


+68.87 


-270 


37.7 


1 


J1245+3356 


133.72 


+83.06 




48.2 


1 


J1151+5437 


140.61 


+60.39 


-155, -126 : 


12.0 


3 


J1208+4540 


144.63 


+69.62 




14.4 


3 


J1220+3853 


149.71 


+76.59 




61.5 


1 


JlOOl+5944 


152.44 


+46.53 


-128 


31.1 


2 


J0928+6025 


154.10 


+42.44 


-141 


70.6 


1 


J1211+3657 


161.27 


+76.95 




36.0 


2 


J0226+0015 


166.57 


-54.38 




37.1 


1 


J0950+4831 


168.96 


+49.10 


-116, +92 


40.0 


1 


J1016+4706 


169.03 


+53.74 


-100 


43.9 


1 


J0212-0737 


171.06 


-62.64 


-91 : 


29.2 


2 


J1103+4141 


172.52 


+63.52 


-120 


44.2 


2 


J0929+4644 


172.65 


+46.02 


-100 : 


25.9 


2 


J0809+4619 


173.33 


+32.21 




31.4 


2 


J0925+4535 


174.42 


+45.54 


-108 


29.9 


2 


J0235-0402 


174.46 


-56.16 


-100 


23.3 


3 


J0803+4332 


176.40 


+31.00 


-95"= 


57.2 


1 


J0843+4117 


180.21 


+38.06 


-114 


59.2 


2 


J1210+3157 


181.58 


+79.96 




55.9 


2 


J0949+3902 


183.55 


+50.57 


-120 


37.9 


2 


J1112+3539 


184.71 


+67.34 


-99"= 


59.5 


1 


J0751+2919 


191.34 


+25.05 


+175 


13.5 


3 


J1104+3141 


195.45 


+66.24 




41.1 


2 


J0912+2957 


195.86 


+42.34 


-100 : 


71.5 


2 


J0401-0540 


196.42 


-39.97 


-99 : 


52.8 


1 


J0914+2823 


198.13 


+42.45 


-110 :,+95 : 


47.2 


1 


J0820+2334 


199.80 


+29.42 




42.0 


1 


J1204+2754 


206.11 


+79.62 




45.4 


2 


J1117+2634 


209.19 


+69.21 




45.9 


2 


J1059+2517 


210.82 


+64.96 


+ 102^^ 


46.0 


2 


J1207+2624 


214.56 


+80.11 




51.4 


2 


J0826+0742 


216.89 


+24.66 


+125 


46.5 


2 


J0910+1014 


219.82 


+35.48 




94.0 


1 


J0947+1005 


225.59 


+43.59 




36.4 


2 


J0943+0531 


230.04 


+40.41 


+90 


55.9 


1 


J0935+0204 


232.39 


+36.84 


+ 132 


55.4 


1 


J1059+1441 


232.72 


+61.19 




37.2 


2 


J1022+0132 


242.16 


+46.07 


+ 135 


55.3 


1 


J1051-0051 


252.24 


+49.88 


+ 131, +210 


37.8 


2 


J1133+0327 


261.35 


+59.87 


+ 102 


52.3 


1 


J1157-0022 


275.71 


+59.64 




49.2 


1 


J1233-0031 


293.11 


+61.99 


+ 130 :" 


29.0 


1 


J1342-0053 


328.82 


+59.37 


-95 


35.6 


1 


J1342+0505 


333.75 


+64.84 


-103= 


38.0 


2 


J1437-0147 


348.72 


+51.37 


-107 


15.4 


3 



Note: Absence of value in the velocity column indicates that there 
is no evidence of HVC along the line of sight. Uncertainty in 
the velocity is dominated by the COS wavelength cahbrations, 
yielding an error of '■^10 kms~ ; the statistical error is generally 
< 3 km s~ , except for those marked with a colon. A velocity 
value foUov/ed by a colon is uncertain owing to blending with the 
lower velocity gas. 

(a): 3iT equivalent width detection limits. Wc estimated W^^j^ near 
the Si III absorption, but the distribution of M^iim for C II or Si II 
would be very similar (although the actual value of Wiim of C II 
or Si II may be different for a given sightline; see ^2.4|l. (6): HST 
Program: (1) 11598, (2) 12248, (3) 11741. (c) HVC only detected 
in Si III and C IV (not Si II or C II). 



Name 


I 


b 

(°) 


.,HVC 

"LSR 

(kms-1) 


1ES1553+113 


21.91 


+43.96 




PG1259+593 
SBS1122+594 


120.56 
141.80 


+58.05 
+54.71 


-142 
-95 


VIIZW244 


136.66 


+32.67 


-93 


3C263 


134.16 


+49.74 


-166, -122 


PG1115+407 


172.23 


+66.67 


-110 : 


MRK421 


179.83 


+65.03 


-119 


HS1102+3441 
TON580 


188.56 
194.94 


+66.22 
+72.03 




HE0238-1904 


200.48 


-63.63 




PKS0405-123 


204.93 


-41.76 


+ 110 


HE0226-4110 


253.94 


-65.78 


+ 156, +208 


HB0435-5304 


261.02 


-41.38 


+ 100 : 


HB0439-5254 


260.69 


-40.90 


+ 106, +300 


RXJ0439. 6-5311 


261.22 


-40.93 


+ 116, +300 


HB0153-4520 


271.80 


-67.98 


+ 110, +197 


PKS2005-489 


350.37 


-32.60 




RXJ2154. 1-4414 


355.18 


-50.87 





Note: A velocity value followed by a colon is un- 
certain owing to blending witli the lower velocity 
gas. 

stellar and AGN TT samples can be directly compared as 
we observe no HVC with W\ < (Wiim) + 2ad. 

Several AGN in CSG09 have spectra with S/N similar 
to the average S/N in the TT sample, but several of them 
have also higher S/N (see Figs. 1 and 2 in CSG09). In Fig.|4] 
we show the measured Si III equivalent widths by CSG09 in 
the individual components of the HVCs of their sample. Fol- 
lowing our guidelines, we removed any HVC that is identified 
only with Si III. Except for 3 high-velocity components, all 
the detected HVCs in CSG09 have also an equivalent width 
larger than (VKiim)- Another way to consider the CSG09 sam- 
ple is for each sightline to retain only the strongest HVC ab- 
sorption components. This is shown by the black histogram 
in Fig. H] For either samples, we find that most of the de- 
tected HVC absorption features have Wx > (Wiim) +2a and 
is larger than {Wnm). 

The sensitivity of the GTO sample is typically similar to 
the best S/N in the TT sample (i.e., comparable to data in 
the GO program 11741). Yet, we estimate equivalent widths 
of the HVC absorption that are similar to those derived in 
the TT sample. Therefore despite a non-uniformity in the 
S/N, the TT, GTO, and CSG09 samples can be merged and 
can be compared to the LHll stellar sample. Fig. H] shows 
that the detections are generally of high significance as oth- 
erwise weak absorption of Si III, Si II, C IV, C II would 
have been detected more frequently in the sample with bet- 
ter sensitivity. This implies that low total H column HVCs 
are not important and that the present sample is essentially 
complete. 



3 COVERING FACTOR AND DISTANCE OF 
THE HVCs 

The sky covering factor of HVCs is a crucial quantity, as it 
is required to statistically determine the distance of the UV 
absorption-line selected population of HVCs (LHll) and for 
estimating the total mass of the HVCs (M ex fcd^, where 
fc is covering factor). Following LHll, we define two classes 
of HVCs: the HVCs with 90 sC |ulsr| < 170 kms"^ and the 
VHVCs with |«lsr| > 170 kms~^. The reason for this sepa- 
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Table 3. Covering factors for the HVCs and VHVCs 



100 200 300 

W^(mA) 

Figure 4. (a) Distribution of the HVC equivalent widths and 
equivalent width limits. The blue filled histogram represents the 
3cr equivalent width detection limit of the AGN TT sample, with 
the vertical red dashed line being the median (42 mA) of this dis- 
tribution. The red hashed histogram is the C II A1334 equivalent 
width of the HVC detected in the TT AGN sample. The black 
hashed histogram is the C II A1334 equivalent width of the HVC 
detected in the LHll stellar sample, (b) Same as (a), but the red 
hashed histogram is the Si III A1206 equivalent width of the HVC 
detected in the CSG09 AGN sample and black hashed histogram 
shows the distribution of the strongest HVC component along 
each line of sight of the CSG09 AGN sample. 



■m /,(HVC) /c(VHVC) /,{HVC + VHVC) 

{%) (%) (%) 


TT AGN Sample with Wnm > 10 mA 


allfe 
fe>0° 


67 58 ±6 %t\ 63 ±6 
60 58 ±6 3^2 62 ±6 


TT AGN Sample with 10 ^ VKum =£ 42 mA 


AUfe 
fe>0° 


34 59 ±9 6l^ 62 ±9 
31 58 ±9 7l^ 61 ±9 


TT+GTO+CSG09 AGN sample 


allfe 

fe>0° 

fe<0° 


133 59 ±4 20 ±4 64 ± 4 
105 57 ±5 13 ±3 61 ±4 
28 68 ± 9 46 ± 10 75 1^ 


LHll stellar sample 


all fe, d > 3 kpc 28 50 ± 9 +g 50 ± 9 
all fe, d > 5 kpc 20 60 ± 10 to 60 ± 10 



Note: m is the size of the considered sample. Errors are 68.3% 
level confidence interval estimated from the Wilson score. 
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ration in velocity is that HVCs were detected toward stars, 
but no VHVC absorption was found in the stellar spectra; 
the precise velocity cutoff is so mew hat arbitrary (as is the 
lower bound at 90 kms~^, see S 3.4 'or more details). In Ta- 
ble |3J we summarize the values of /c for different AGN sam- 
ples and the LHll stellar sample (see LHll and Q. In this 
table, m is the sample size (i.e., the number of sightlines). 
We assess a 68% confidence interval for each value of /c 
using the Wilson score interval for a binomial distribution. 
The main results can be visualized in Fig. |5] are listed in 
Table [3] and can be summarized as follows: 

(i) /c does not strongly depend on the sensitivity of the 
present data; 

(ii) Comparing the H I surveys with ours, we conclude 
that 74% of the HVC-hVHVC directions have iV(H I) < 
3 X 10^** cm"^ and 46% have iV(H I) < 8 x lO" cm"^ (panel 
(a) in Fig. §; 

(iii) the HVC detection rates in the stellar and AGN sam- 
ples overlap within Icr, strengthening the LHll conclusion 
that most of the HVCs are near the Milky Way disk, within 
about heliocentric distances 5-15 kpc (panel (b) in Fig. pi; 



Figure 5. Covering factors for the HVCs and VHVCs: (a) com- 
parison between the H I emission (A'^(H I) limited) and UV ab- 
sorption (any Ar(H I), i.e., NiJA I) < IQis cm'^ to Ar(H I) > 
10^" cm~^) samples; (b) comparison the covering factors de- 
termined from the TT-I-GTO-I-CSG09 AGN and LHll stellar 
(where two different heliocentric distances cutoff are considered 
at 3 and 5 kpc) samples; (c) VHVC covering factor from the 
TT-I-GTO-I-CSG09 AGN and LHll stellar samples. 

(iv) the HVC detection rate in the stellar sample drops 
with decreasing distances (panel (b) in Fig. [5|; 

(v) VHVCs are far more frequent at 6 < 0° than at fe > 
0° , which is largely due to large Galactic sky covering of the 
Magellanic Stream at fe < 0° (panel (c) in Fig. |5|. 

In the following sections, we explain and discuss these con- 
clusions. 



3.1 HVC covering factor from the AGN samples 

We first consider solely the TT sample because we are cer- 
tain that no bias was introduced in selecting the targets that 
would favor or disfavor the detection of HVCs. In order to 
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confirm that our HVC samples are essentially complete, we 
consider one sample with VKum 1? 10 mA (i.e., the entire TT 
sample) and a sample with 10 ^ VKum ^ 42 mA (i.e, a sam- 
ple with a higher sensitivity). The results are summarized 
in Table [3] The entire and more sensitive TT samples show 
that /c is essentially independent of the equivalent width 
sensitivity of the present data (a sample with Wum ^ 42 mA 
gives /c = 56%, similar to the /c-values listed in Table ^. 
As demonstrated in Fig. |4] and discussed in §2.41 there is 
no evidence for a large population of weak metal absorbers 
with a significant fc- 

Fig. [1] shows that there is no HVC at fe > 75° in the TT 
sample, which is also confirmed with the combined sample 
(see Fig. [6|. For 30° < 6 < 75° the HVC detection rate is 
necessarily higher at 69% compared with 58% for the 6 > 0° 
sky. The absence of HVCs at fe > 75° could be due in part 
to statistical fluctuations over the small solid angle covered 
by these latitudes. If the HVCs follow a thickened disk-like 
distribution, a deficit of clouds at fe > 75° would also be 
more likely. 

In Fig. [6] we show the distribution of sightlines and 
HVC detections for the combined sample that includes the 
TT, GTO, and CSG09 samples. This combined sample has 
much better coverage at fe < 0° where it is larger than the 
TT sample by a factor 4. The observations are still heav- 
ily weighted toward the northern Galactic sky. In Table |3J 
we summarize the detection rates of the combined sample. 
We again consider the HVC, VHVC, and HVC+VHVC cat- 
egories in 3 sub-samples: the entire Galactic sky, the north- 
ern Galactic hemisphere (fe > 0°), and the southern Galactic 
hemisphere (fe < 0°). 

For b > 0°, the hit rates are very similar for the HVC 
and HVC-I-VHVC between the TT and full samples. Within 
the errors the VHVC hit rate is also similar and small at 
fe > 0° . Comparing the entire sky yields similar conclusions, 
but again the sample is heavily weighted to fo > 0° with 
~ 4 times more AGN than at fe < 0°. Given this similarity, 
we will adopt the results from the combined sample for the 
remainder of this work. There is, however, a clear difference 
at fe < 0°: the covering factor of the VHVC is substan- 
tially larger, with a 46% covering factor compared to 13% 
at fe > 0° (see Table [3] and Fig. [5|. There is also a tentative 
increase in /c(HVC) at fe < 0° compared to that at 6 > 0°, 
although the values overlap with the errors. As we discuss 
further in §3.4| the southern Galactic sky is dominated by 
the Magellanic Stream ( [Putman et al.pQQ S: Nidc ver et al.| 
20081 ), and this mostly explains the considerable difference 
in /c(VHVC) between the two Galactic hemispheres. We 
defer a discussion on the sky coverage of the VHVCs and 
implications for the origins and distances of the VHVCs to 

The covering factors derived from our samples using 
O I, C II, C IV, Si II, Si III, Si IV etc., are quite consis- 
tent with those determined for HVCs with O VI absorption, 
which are in the range of 60%-75% ( [Sembach ct al. 2003 ; 
Fox et al.|2006 l. The fact that the covering factors are simi- 
lar is not surprising in view of the multiphase nature of the 
HVCs, which often show absorption by O I, C II, C III, C IV, 
O VI at similar velocities in a given sightline (e.g., Ganguly 



0° (B *^0 8 O 



-^ 



Figure 6. Aitoff projection of the combined (TT-(-GTO-(-CSG09) 
sample. A filled circle means at least one HVC with |^'LSR| > 90 
kms~^ is detected along the line of sight. An open circle means 
that no HVC is detected along the sightline. 



face between the hot Galactic corona and the cooler ionized 
and neutral HVCs (e.g., [Sembach et al.||200^ [Collins et al.| 
2005}|Kwaket al.|2011| ). 

We note that the difference in the hit rates between our 



et al.|2005| [Collins et al.|2005||Fox et al.|2004|r2006 , Zoch et. 



al.[[2008 L As warm HVCs move in the Galactic halo, some 
of the "O VI HVCs" are likely the coUisionally ionized inter- 



analysis and CSG09 (and [Shull et al. 2009 where they found 
/c(HVC + VHVC) = 80-90% based on the search of Si III 
HVC absorption in AGN spectra) is due to three causes: i) 
there was a bias in the CSG09 sample where AGN initially 
targeted for HVC science were not removed, which artifi- 
cially increased the covering factor; ii) while CGS09 con- 
firmed in many cases the HVC absorption with ions other 
than Si HI, in several cases they rely solely on the Si III 
absorption and assumed that the intergalactic contamina- 
tion is small; and iii) as we have just demonstrated, the TT 
and combined samples are not uniformly distributed in both 
hemispheres, thus our sample weights the northern hemi- 
sphere more strongly than CSG09 (and the CSG09 sample 
is more balanced between the two Galactic hemispheres) . A 
uniform measure of the entire Galactic sky covering factor 
of the HVC-(-VHVC would likely be somewhat higher than 
the 64% derived here. In view of this imbalance, a straight 
average between the two Galactic hemispheres is more ap- 
propriate, yielding a covering factor for the entire sky of 
/c(HVC + VHVC) = 68% ± 4%. That is the value adopted 
in Fig. [5] (panel (a)). 



3.2 Distance of the HVC population 

Having set a robust HVC covering factor of /c(HVC) = 
59% ± 4% toward the AGN sightlines, we can revisit the 
LHll method to statistically determine the distance of the 
HVC population by comparing the HVC covering factors 
determined toward sightlines with no distance limits (AGN) 
and with distance limits (stars). With their stellar sample, 
LHll found /c(HVC) = 50% ± 9% (m = 28). The /^-values 
of the AGN and LHll stellar samples are therefore consis- 
tent within the 68% confidence interval (see also Tablelsland 
Fig. [5|. The same ions and techniques were used to search 
for the HVC absorption in the stellar and AGN spectra, al- 
though finding HVC absorption in the stellar sample often 
relied only on C II and Si II (Si III was often too noisy and 
C IV could be blended with the stellar absorption). While it 
may be possible for more highly ionized HVCs to be missed 
in the stellar sample, this is unlikely to be important since 
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a large majority of the HVCs would still show absorption in 
C II and Si II as discussed in |2.3[ 

Although the results are consistent within a 68% confi- 
dence interval, the difference between the stellar and AGN 
fc values is large enough that it could suggest a fraction of 
HVCs being at much greater distances if fc is not distance- 
dependent. However, owing to interstellar drag on these 
clouds, HVCs should slow down as they fall onto the Milky 
Way disk ( Benjamin fc Danly||1997 l. More recent hydrody- 
namical simulations have also followed the infall of HVCs 
through the Galaxy and predict that the HVCs become dis- 
rupted by their interaction with the corona (e.g., |Bland- 



Table 4. Stellar sample with stars at d > 2 kpc, \z\ ^ 1 kpc 



Name 



Hawthorn et al.||2007| [Peek et al.|[2008| [Heitsch fc Putman 



2009[|Kwak et ai.|2011[|Joung et al.|2012[ ). These two argu 
ments imply that the number of HVCs should decrease at 
smaller d or z. As LHll selected their stellar sample so that 
the stars are at \z\ > 3 and d > 3 kpc (only one star was 
at z = —2.6 kpc, but d = 8.7 kpc), we searched MAST for 
stars at smaller distances d > 2 kpc (setting this limit to 
avoid probing only our local environment) and [zj > 1 kpc 
(and |6| > 15° to avoid probing HVCs linked to phenomena 
occurring in the Galactic disk, see, e.g., Lehner et al.||2011 l 
to study the dependence between fc and d. We found an 
additional 8 stars that satisfy these criteria with the appro- 
priate coverage of at least C II A1334, Si II A1260, Si III 
A1206. We summarized the full (LHll plus 8 stars) stellar 
sample in Table |4] where the stars are ordered in z-height 
from the Galactic plane. 

With this sample of stars, we estimate the values of fc 
as function of d or \z\, i.e., we calculated fc for stars with 
d > 2 kpc (at any |z;| > 1 kpc), > 3 kpc, etc., and similarly 
with l^l (and any d > 2 kpc). The results are summarized 
in Fig. It] where there is a distinct trend between fc and d, z. 
At d < 3-4 kpc {\z\ < 2-3 kpc), fc is smaller and not consis- 
tent within the 68% confidence interval with fc determined 
from the AGN sample, but as d (|2:|) increases, fc increases, 
and plateaus at d > 4 kpc {\z\ > 3 kpc) to the fc value 
determined from the AGN sample. For d Js 5 kpc (|2| > 3-4 
kpc), fc = 61% ± 10%, in remarkable agreement with the 
result from the AGN sample. Based on this result, if we set 
a cutoff at z = 4 kpc, we find: 

- |0| ^ 4 kpc, fc = 58% ± 11% ((d) = 11.5 ± 6.5 kpc); 

- |«| < 4 kpc, fc = 20%+^^f° ((d) = 4.1 ± 1.8 kpc). 

A larger sample of stars would be needed to better char- 
acterize the infall of HVCs onto the Galactic disk, but our 
results strongly suggest that the population of HVCs de- 
creases with decreasing \z\ (and d), and hence the HVCs may 
either be disrupted or decelerated to become intermediate- 
velocity clouds (IVCs, 50 ^ fLSR < 90 kms"'^) or low- 
velocity clouds. A detail description of the (UV-selected) 
IVC population is beyond the scope of this work, but it 
would be extremely valuable to undertake to understand 
how HVCs may reach the Galactic disk. 



The excellent agreement in the HVC detection rates be- 
tween the stellar and AGN samples strengthens the conclu- 
sions reached by LHll: most of the HVCs are within 5-15 
kpc of the sun. The drop in the population of HVCs with \z\ 
is consistent with models predicting their fall through the 
Galactic halo onto the disk. All this implies that the HVCs 
are distant enough to have sufficient mass and close enough 



I b d \z\ ^HVC 

n n (kpc) (kpc) (kms-i) 



HD195455'' 


61 


-27 


2.2 


1.0 


+94 


HD116852'^ 


304 


-16 


4.5 


1.3 




HD215733'" 


85 


-36 


2.9 


1.7 


... b 


HD149881'' 


31 


-1-36 


2.9 


1.7 




JL212'' 


303 


-61 


2.2 


1.9 




NGC6723-III60 





-17 


8.7 


2.6 


-90 


HD100340'" 


258 


61 


3.0 


2.6 




HD18100" 


217 


-62 


3.1 


2.7 




HD119608'' 


320 


43 


4.1 


2.8 




HD121968 


334 


-1-55 


3.8 


3.1 




PG1511-I-367 


59 


-t-59 


3.8 


3.2 




NGC104-UIT14 


306 


-45 


4.5 


3.2 




HD233622 


168 


4-44 


4.7 


3.3 




EC10500-1358 


264 


+40 


5.2 


3.3 


+97 


PG1704-I-222 


43 


-f32 


6.9 


3.6 




PG0855-I-294 


196 


+39 


6.5 


4.1 


+93, +107 


PG2219-I-094 


73 


-40 


6.6 


4.2 




PG0832-I-675 


148 


+35 


7.5 


4.3 


-123 


PG0955-I-291 


200 


+52 


5.5 


4.3 




NGC6341-326 


68 


+35 


8.2 


4.7 


-94; 


NGC6205-Barnard29 


59 


+41 


7.7 


5.0 


-107 


NGC5904-ZNG1 


4 


+47 


7.2 


5.3 


-140,-120 


PG1610-I-239 


41 


+45 


8.4 


5.9 




HS1914-I-7139 


103 


+24 


14.9 


6.0 


-175,-118 


PG0122-I-214 


133 


-41 


9.6 


6.2 


-160,-91 


PHL346 


41 


-58 


8.7 


7.4 




SB357 


301 


-81 


7.9 


7.8 




PG0914-I-001 


232 


+32 


16.0 


8.4 


+100, +170 


PG0009-I-036 


105 


-58 


10.8 


9.1 




vZ1128 


43 


+79 


10.2 


10.0 




PG1708-I-142 


35 


+29 


21.0 


10.0 




PG1002-I-506 


165 


+51 


13.9 


10.8 


-102, +101 


NGC5824-ZNG1 


333 


+22 


32.0 


12.0 


-160: 


PG1323-086 


317 


+53 


15.8 


12.6 


-91 



Note: a: Stars not originally in the L Hll stellar sample, b: 
There is an IVG at iilsr = —84 kms~^ [Fitzpatrick &: Spitzer| 
119971 A colon means the result is tentative in view of the 
low signal-to-noise of the data (removing these stars from the 
sample would not change our conclusions). 



that their infall rate balances the star formation rate in the 
disk of the Milky Way (see LHll and ^5^. In plwe also 



show that several of the ionized HVCs are closely related to 
the large H I complexes given their proximity on the celestial 
sphere and in three dimensional space. The (neutral and ion- 
ized) HVCs are therefore very likely a large part (if not all) 
of the source of gas required for continued star formation. 



3.3 HVC covering factors from H I emission and 
UV absorption data 

Searching for high- velocity H I emission with i?;LSRJ J? 90 
kms~^ using the LAB survey toward the sightlines in the 
combined sample, we find fc — 18% for logA^(H I) ^ 18.5. 
Given that more sightlines in the CSG09 and GTO samples 
go through large H I HVC complexes (Magellanic Stream, 
complex C and complex A), it is not surprising that the H I 
LAB HVC detection is higher in these samples {fc ~ 22%) 
than in the TT sample {fc ~ 12%). The detection rate from 
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niques are different and the absorption-line method samples 
comparatively sparsely the Galactic sky. Our AGN sample 
and the absorption-line technique are well suited for finding 
HVCs with both high and low H I column densities. 

Deeper radio H I emission surveys reached a sensi- 
tivity of logAf(H I) ~ f7.9 (, Murphy et aL||1995| [Lock 
man et al. 20021. At this limit, 37% of the sky is cov- 



12 3 4 5 

zl (kpc) 

Figure 7. The HVC covering factors toward the stars in the 
stellar sample summarized in Table |4] which includes the LHll 
stellar sample plus 8 stars at 1 < 1^1 < 3 kpc. Top panel: fc vs. 
d. Each data point corresponds to the value of fc for stars at 
distances > di, where di = 2, 3, ... kpc. The horizontal dotted line 
is the HVC covering factor determined toward the AGN. Bottom 
panel: same as top panel but fc is shown against z. As stars with 
smaller d or \z\ are removed from the sample, fc converges on the 
/c-value determined from the AGN sample. 



the combined sample is entirely consistent with that derived 
for the entire radio H I sky survey at a similar sensitiv- 
ity (Wakker 19911. We emphasize this result as the tech- 



ered by H I HVCs+VHVCs, which is still substantially 
smaller than fc in our UV absorption survey. Extrapolat- 
ing to logiV(H I) ~ 17, |Lockman et al.| (20021) argued that 



the HVC+VHVC covering factor would be about 60%, close 
to our findings {fc = 68%). However, they assume that the 
constant covering factor per decade of logA''(H I) observed 
between 19 and 17.9 dex applies down to 17 dex. While di- 
rect estimates of A'^(H I) via H I absorption is difficult. Fox 
et al. ( 2006 1 searched for H I Lyman series absorption as- 



est column UV-selected HVCs have logAf(H I) < 14.7-15.0 
while the highest have log A''(H I) > 20; thus there is at least 
a 5 dex range in A'^(H I). The detection of O VI or C III or 
any other metal line species in low H I column density HVCs 
still imphes N{R) > 10"-10^*(Z/Zq)-^ cm'^ Thus, even 
though the H I column density in these HVCs may be small, 
the total hydrogen column densities can still be substantial. 
Comparing the H I surveys with ours, the relative 
covering factors imply about 74% of the HVC-I-VHVC di- 
rections have N(H I) < 3 x 10^* cm"^ and 46% have 
N{}i I) < 8 X 10^^ cm^^. We emphasize again the popu- 
lation of HVCs with extremely low total H column density 



is unimportant (see f2.4| 



3.4 



Covering factor and distribution of the 
VHVCs 



LHll did not find any HVCs with |wlsr| > 170 kms"\ the 
so-called VHVCs, in absorption in the stellar spectra. They 
concluded that the VHVCs must be at distances d > 10-20 
kpc. Our analysis of the covering factors of VHVCs using the 
much larger sample of AGN demonstrates that at b > 0°, 
the covering factor is quite small, /c(VHVC) ~ 13%(see Ta- 
ble [s] and Fig. [5|. In fact the TT sample nearly completely 
missed the VHVCs. The LHll stellar sample and TT AGN 
sample are distributed similarly across the Galactic sky (see 
Fig. [T|, but the LHll sample is about 2.5 times smaller. It 
may not be therefore surprising that no VHVC absorption 
was observed in the spectrum of a star (/c(VHVC) < 13% 
for a 95% confidence interval for the stellar sample). Some of 
the stars cover the fe < 0° hemisphere where the /c(VHVC) 
is much larger, but, as we show below, most of the VHVC 
absorption is associated with the Magellanic Stream. The 
absence of Magellanic Stream absorption toward stars con- 
firms that the Stream is indeed far away. Using our large 
sample of AGN, we can learn more about the origins of these 
clouds as we now demonstrate. 

In Fig. Is] we show the sky distribution of VHVC absorp- 
tion from our combined sample as well as H I 21-cm emission 
from the NRAO 140-f t (logAf(H I) ^ 17.9) observations of 
Lockman et al. (20021. We do not show the lower velocity 



sociated with O VI and/or C III HVCs. Their results show 
that the H I column densities vary tremendously: the low- 



HVCs in order to highlight solely the VHVCs. In this figure, 
open symbols are used if no VHVC is seen in absorption or 
emission. In the bottom panel of this figure, the VHVCs with 
|wlsr,| > 170 kms~^ are shown, while the top shows only the 
VHVCs with |wlsr| > 210 kms~^. The overall sky coverage 
is excellent with only a lack of H I data in the 6 < 0° and 
240° < I < 360° quadrant because this region of the sky 
cannot be observed from Green Bank. The velocity distri- 
bution seen in emission and absorption is strikingly similar, 
which reinforces that the VHVCs seen in emission and ab- 
sorption must be similar structures with different ionization 
levels (see Q. It is therefore not surprising that there is 
a remarkable correlation between these maps and all sky 
H I 21-cm HVC map (e.g., [Wakker et al.||2003[ ), with many 
of the VHVCs superimposed onto or near large known H I 
HVC complexes; the Magellanic Stream at similar velocities, 
complexes C and A, the anti-center HVCs at higher veloc- 
ities. Wakker (20011 and Lockman et al. (20021 attributed 

170 



the gas moving at velocities with - 
kms"^ in the b > 0° and 0° < I 



-210 < 

< 180° 



llLSR < 

quadrant pri- 
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Figure 8. Aitoff projection maps showing the sightlines where a VHVC is detected toward an AGN (filled circles from UV absorption 
TT+GTO+CSG09 data; filled square from the H I 21-cm emission data from Lockman et al. 2002). On the top panel, detections are 
VHVCs with |fLSR.| > 210 kms~ , while on the bottom panel, detections are VHVCs with I^lsrI > 170 kms~^. The region defined by 
6 < 0° and 240° < 6 < 360° cannot be observed from Green Bank, explaining the lack of H I emission data in this region. All the open 
symbols indicate that no VHVC absorption or emission with [tiLSRl > 210 kms~^ {top) or |i>lsrI > 170 kms"'^ [bottom) is found along 
the sightline. 



marily to complex C. In the top panel of Fig. [S] where we 
only show now the VHVCs with JiilsrI > 210 kms~^, the 
VHVCs affiliated with complex C have almost completely 
disappeared (as well as some of the anti-center complex). 
The VHVCs that remain present correlate extremely well 
with the H I map of the Magellanic stream and its lead- 



that the most sensitive H I emission and UV absorption 
data suggest it could be wider (i.e., the Stream has a large 
ionized envelope) and could reach higher positive galactic 
latitudes than previously thought. Some of the VHVCs near 
I = 120°-130° and b = -20° to -30° are likely associated 
with the M31 and M33 galaxies (Thilker et al.||2004[ |West-| 



ing arm (Wakker et al. 2003, Nidever et al. 20081, except meier et al.|2005 Putman et al.|2d09l, but their solid angle 
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Table 5. Heliocentric distance (rfHVc)i height (|zhvcI 
the HVCs and association with H I HVC complex 



■ <^HVcl sin6|), and Galactocentric distance" (Rg) of 



Name 



I 

n 



n 



dnvc 
(kpc) 



IZHVCI 
(kpc) 



Rg 

(kpc) 



"LSR 
(kms-1) 



Associated H I 
HVC Complex 



Stellar Distance 
References 



NGC6723-ni60 





-17 


<8.7 


<2.6 


<0.2 


-90 


None 


(1) 


NGC5904-ZNG1 


4 


+47 


< 7.5 


< 5.5 


< 3.4 


-140 


Complex L? 


(1) 


NGC5904-ZNG1 


4 


+47 


< 7.5 


< 5.5 


< 3.4 


-120 


Complex L? 


(1) 


M13-Barnard29 


59 


+41 


< 7.1 


<4.7 


< 7.4 


-121 


Complex C 


(1) 


HD195455'' 


61 


-27 


< 2.2 


< 1.0 


< 7.7 


+94 


Complex gp? 


(2) 


HS1914+7139 


103 


+24 


< 14.9 


<6.0 


< 17.6 


-118 


Outer Arm/Complex C? 


(3) 


HS1914+7139 


103 


+24 


< 14.9 


<6.0 


< 17.6 


-175 


Complex C 


(3) 


PG0122+214 


133 


-41 


<9.6 


<6.2 


< 14.5 


-91 


Cohen Stream 


(3) 


PG0122+214 


133 


-41 


<9.6 


<6.2 


< 14.5 


-160 


WW503 


(3) 


PG0832+675 


148 


+35 


> 8.1 


^ 4.7 


^ 14.6 


-145 


Complex A 


(4) 


PG0832+675 


148 


+35 


<8.1 


<4.7 


< 14.6 


-123 


Ionized complex A? 


(4) 


PG1002+506 


165 


+51 


< 13.9 


< 10.8 


< 17.1 


-102 


Complex M 


(5) 


PG1002+506 


165 


+51 


< 13.9 


< 10.8 


< 17.1 


+101 


None 


(5) 


PG0855+294 


196 


+39 


<6.5 


<4.1 


< 13.4 


+93 


Cloud WW113 


(6) 


PG0855+294 


196 


+39 


<6.5 


<4.1 


< 13.4 


+107 


Cloud WW113 


(6) 


PG0914+001 


232 


+32 


< 16.0 


<8.4 


< 20.0 


+100 


Complex WB 


(6) 


PG0914+001 


232 


+32 


< 16.0 


<8.4 


< 20.0 


+170 


Complex WA 


(6) 


EC10500-1358 


264 


+40 


< 5.2 


< 3.3 


<9.8 


+97 


Cloud WW95 


(2) 


NGC104-UIT14 


306 


-45 


>4.5 


> 3.2 


> 7.1 


C 


Complex WE 


(1) 


PG1323-086 


317 


+53 


< 15.8 


< 12.6 


<6.7 


-91 


None 


(7) 



(a): Galactocentric distance of the HVCs, where Rg = {Rq + (i^HVC cosb)^ — 2diiYC Rq cosl cos 6)"'^, 
assuming the distance from the Sun to the Galactic center Rq = 8.5 kpc. (b): Not in the LHll sample, (c): 
No absorption is detected, but strong H I emission is observed along this line of sight. References for the 
heliocentric distances of stars that were used to estimate the distance limits of the HVCs: (1) Harris 1996] 
2010| (2) [RoUeston et aL]|1997| (3) [Ramspeck et al.||2001| (4) |Wakker|[200l1 (5) [Ringwald et al.jil998| (6) 
RoUeston et al.|1999| (7) |Moehler &: Heber|1998l 



is extremely small. We therefore conclude that most of the 
VHVCs with IulsrJ > 210 kms^'^ are dominated by the 
Magellanic Stream at 6 < 0° and its leading arm at 6 > 0°. 
Only a very few of the VHVCs, such as those associated 
with M31 and M33, may not be part of the Stream. Inde- 



pendently, Putman et al. (2011 1 showed recently that many 



compact VHVCs have a head-tail structure (a signature of 
interaction with a hot diffuse medium) and can be associ- 
ated with the Magellanic Stream. From our discussion, these 
small head-tail H I VHVCs must have larger envelopes of 
photoionized gas and collisionally ionized gas. The VHVCs 
are therefore not randomly distributed on the Galactic sky 
and appear mostly to be associated to known large struc- 
tures. However, while our AGN sample is large, the Galactic 
sky can only probed sparsely along pencil-beam sightlines 
with the absorption technique. It is therefore quite plausi- 
ble that we might miss a population of compact HVCs with 
very small angular scales at very large distances (or that 
they may be confused with the more nearby VHVCs) as 



those possibly found with the ALFALFA survey ( Giovanelli 
et al.|[2010l. The compact HVCs are therefore the only re- 



fying their basic physical properties. The distances to many 
H I HVCs have been determined using stars in direction of 
the H I contours with N{Ri) ^ 3 x 10^* cm"^. The high- 
velocity gas seen in absorption toward AGN outside the H I 
contours has previously been associated with the extended 
HVCs seen in 21-cm emission based on their angular proxim- 
ity and similarity in velocities (e.g., [Tbx et al.|["2004| |Collins| 
et al.|2005|. With the LHll stellar sample, we can go a step 



maining candidates for dark matter mini-halos in the HVC 
population. 



further and directly compare the distances of H I HVCs and 
HVCs seen in UV absorption. We emphasize that for most 
of the LHll stellar sightlines, no H I 21-cm emission (at the 
LAB survey sensitivity) is observed. 

To search for H I HVCs associated with the UV- 
selected HVCs, we used the H I maps in |Wakker et al.| 
( 2003 1 and Wakker ( 2001 1 as well as information provided by 
B.P. Wakker (private communication, 2011). In Appendix [A] 
we summarize in detail the relation between the HVCs seen 
in UV absorption for each sightline with the H I com- 
plexes, while in Table [5J we summarize the HVC distance 
constraints from the LHll stellar sample for which the de- 
tection is unambiguously a high- velocity interstellar feature 
and the S/N was high enough for their reliable detectionrl 
In this table, we also added the additional star at smaller 



4 ASSOCIATION OF THE NEUTRAL AND 
IONIZED HVCs IN 3D SPACE 

Determining the distance of the HVCs has been critical for 
associating HVCs with flows occurring near the Milky Way 
rather than in the IGM of the Local Group and for quanti- 



■^ Specifically the following 4 stars were not considered: 
NGC6341-326 and NGC5824-ZNG1 owing to the low S/N in their 
COS spectra, and PG1243+275 and PG0934+145 for stellar con- 
tamination as noted by LHll. For the latter star, there is a possi- 
bility that both stellar and interstellar high-velocity components 
are present at the same velocity. 
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Figure 9. Normalized profiles of interstellar C II and Si II absorption seen in tlie spectra of several stars. The HVC components are 
marked by the vertical dotted lines. Other features are either stellar or interstellar. The spectra were normalized within ±600 kms~^ 
from the absorption lines using Legendre polynomials with a degree dp (see, e.g., |Lehner et ar||201l| . For many stars, a fit with a low 
degree, dp = 1-3, was used, but some stars (e.g., PG0855+294, EC10500-1358) have complicated stellar profiles, especially near C II, 
requiring higher values of dp > 6. There is a broad stellar line next the interstellar C II feature at positive velocities in the spectrum of 
EC10500-1358 that is only approximately modeled (a feature that is similar to that seen in PG0122+214 that we did not attempt to 
remove in this case as no HVC is seen at positive velocities in this direction). 



d with a HVC detection (see ^3^ and Table |4|. The he- 
liocentric distance, z-height above and below the Galactic 
disk, Galactocentric distance, velocity, and possible associa- 
tion with H I HVC complexes are listed in Table [5] We show 
the normalized interstellar absorption spectra in Figs. |9] |10[ 
fTTlfor those stars whose STIS or COS spectra have not yet 
appeared in the literature (see Appendix [A|). We refer the 
reader to LHll for more information about the data acquisi- 
tion and reduction of the stellar sample. We emphasize that 
many of these HVCs are also detected in other transitions 
and/or species (e.g., O I, Al II, Si III, C IV, see figures in 
Lehner & Howk 2010, LHll, Zech et al. 2008, and Fig.llOl. 



From our detailed discussion in Appendix lA] and the 
results summarized in Table [5] we find that the distances 
to many of the predominantly ionized HVCs seen in the 
stellar spectra are consistent with the distances found for 
the H I HVCs nearby on the sky at similar LSR velocities. 
This implies that many of the HVCs seen in absorption are 
the extended diffuse envelopes of the denser clouds seen in 



H I emission. In view of the much larger covering factor of 
the predominantly ionized HVCs and the proximity of the 
sightlines to many H I complexes, all the H I HVCs are 
likely to have extended ionized envelopes. However, in both 
the stellar and AGN sightlines there are also ionized HVCs 
where no H I complex is observed within 20°-30° of the 
sightline, implying that the ionized high-velocity gas is not 
only found near H I complexes but also in regions devoid of 
H I 21-cm emission (at least for a H I sensitivity of about 
> 3 X 10^* cm-2). 



5 IMPLICATIONS 

5.1 Origins of the gas flows in Milky Way halo 

The results discussed in §3.2| and !|4]have demonstrated that: 

(i) most of the HVCs are within about heliocentric dis- 
tances 5-15 kpc. 
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Figure 10. Normalized absorption profiles of various metal- 
lines and LAB HI emission line profile {top panel) toward 
PG0832-I-675. The HVC components are marked by the verti- 
cal dotted lines. The main H I emission is observed at —145 
kms~ (note that the most negative emission component is not 
observed along this sightline using a 1' beam, Wakker et al. 2003) 
and is part of complex A. The main interstellar absorption is 
seen at —123 kms~^, although some absorption is also seen at 

— 145 kms"'^. The stellar spectrum is complicated with many 
stellar feature (the stellar velocity is —70 kms~^ where absorp- 
tion is seen in Al III as well as Si IV and C IV). We argue that 
PG0832-I-675 must be at the front edge of complex A and the 

— 123 kms~ component could be part of its ionized envelope. 



(ii) the population of HVCs drops with decreasing dis- 
tances or z-height, and 

(iii) several HVCs seen in UV absorption are the extended 
ionized halos of their neutral counterparts. 

The first item implies that HVCs represent flows occurring 
in the Milky Way halo within about one Galactic radius from 
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Figure 11. Normalized profiles of interstellar C II, Si II, Si III 
absorption seen in the STIS E140H spectrum of HD195455, which 
was not originally in the LHll sample. The component at +70 
km s~^ is consistent with the velocity of the IVC complex gp. The 
-f94 kms~^ component might be a higher velocity component of 
this complex. 



the Galactic center. The second item is consistent with in- 
fall of HVCs onto the Milky Way disk. With the third item, 
it is not a big leap to conclude that the metallicities UV- 
selected HVCs are similar to those of their H I counterparts. 
This is critical, as the metallicity can be difficult to estimate 
for the ionized HVCs, often requiring large ionization cor- 



rections (although see Zech et al. 20081. Past and current 



estimates of the metallicities of HVCs find a wide range of 
values from about a few percent solar to super-solar metallic- 
ity (e.g., |Wakker|2001| [Tripp et al. 1 20031 IColUns et al.|2003 



Fox et al.|2004||Zech et al.|2008||Shull et al.|2011||Yao et ah] 
2011| [Tripp fc Song||2012| B.P. Wakker, 2011, private com- 
munication), with a rough average around 0.2-0.4 Zq. As 
for their H I counterparts, several ionized HVC complexes 
must therefore have an extragalactic origin, i.e., not all are 
recycled gas recently ejected from the Milky Way (although 
we note that if some HVCs were ejected from regions situ- 
ated at Galactocentric radii beyond 10-15 kpc, a metallicity 
of 0.4 solar would not be unjustified based on the Galactic 
abundance gradient, e.g., Chiappini et al.[|2001 l. 

While some of these HVCs may be the result of cold 
stream accretion wherein metal-poor gas flows onto galax- 



ies along dense intergalactic filaments (Keres & Hernquist 



2009 Putman 2006 1 or remnants clouds from the Galaxy for- 
mation ( [Oort|1970[[Peek et al.|2008| ), it is also possible that 
much of this gas may come from the interactions with nearby 
dwarf galaxies and their outflows. Such a scenario may be 
seen farther away with the Magellanic Stream and Large 
Magellanic Cloud outflows. The masses of the Stream and 
outflows from the LMC are quite substantial ( Putman|2006 



[Nidever et"aLl[2008[ [Lehner et aTlfSoogj ). Being so far away, 
it is unclear if these HVCs will be able to reach the Galactic 
disk, instead they may mostly enrich the hot Galactic corona 
( [Fox et al.|2010[[Putman et al.|2011[ ). However, it is not in- 
conceivable to imagine that closer dwarf galaxies may feed 
the Milky Way disk via their outflows and tidally stripped 
material. One way to differentiate these scenarios would be 
to study the distribution of the HVC metallicities. At higher 
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redshifts H I absorbers with 16 sj logAf(H I) < 18.5 that 



trace the circumgalactic gas of galaxies (see f 5.3 1 are some- 



times observed at a metallicity of about 2% solar or less 
(Tripp et al.||2005| IRibaudo et al.||2011l N. Lehner et al. 



2012, in preparation), and these may have originated from 
cold flow accretion ( Ribaudo et al.|2011 1. HVCs with metal- 
licity < 2% solar (hke perhaps complex A, see Wakker|2001 1 
could be candidates for cold stream accretion. On the other 
hand, for > 20% solar metallicity (like complex C), the gas 
has been polluted to levels as observed in the LMC outflows 
or Stream, which could be a signature of interaction between 
nearby dwarf galaxies and the Milky Way or of pollution of 
more metal poor gas by fountain or thick disk gas. 

There are also HVCs with solar-like metallicity, sugges 
five of galactic fountain flows (e.g. 



observations. The mass of ionized gas derived from these ob- 
servations is similar to the H I mass within the H I contours 
(i.e., where logA''(H I) > 18.5), but the Ha emission con- 
tours extend far beyond the H I contours (Hill et aL]|2009 



Barger et al.||20i2l, implying a greater mass of ionized gas 



While there is still some uncertainty in the fraction of 
the HVC population that is infalling, the metallicity, and the 
infall timescales, LHll conservatively derived an infall rate 
of 0.4-1.4 Mq yr~^. This is enough to balance the present- 
day total star formation rate (1.9±0.4 Mq yr~^) of the Milky 



Way ( Chomiuk fc Povich|20li| | taken into account that the 
amount of infall needed to sustain star formation may be 
smaller than the star formation rate according to chemical 
evolution models (Chiappini et al. 2001) and stellar mass 



Bregman 1980 Fraternali & Binney 20061. For example, 



Shapiro fc Field||1976{ loss models ( [Leitner fc Kravtsov|2011[ ) 



Zech et al. ( 2008 1 derived a super-solar metallicity for two 



of these HVCs toward the Galactic center. This region is es- 
sentially devoid of H I emission at high-velocity (except for 
complex L), but high- velocity absorption is observed toward 
several AGN and stars ([Sembach et al.|2003||Fox et al.|2006 



IKeeney et al.|2006|[Z"ech et al.|2008||Bowen et al.|2008[ this 
paper, see Appendix|A[| . The high metallicity and absence of 
neutral gas are consistent with the HVCs tracing large-scale 
recycling flows within the central region of the Milky Way 
as those predicted by Galactic fountain models. 



5.2 Fueling the Milky Way 

As the HVCs with 90 ^ |ulsr| < 200 kms"^ are within a 
few kiloparsecs from the Milky Way irrespective of their H I 
content, they are also the most likely source of gas for fu- 
eling continued star formation in the Milky Way. |Collins et] 
al. ( 2005 1 and others argued that the segregation in positive 



and negative radial LSR velocity of HVCs with the Galactic 
coordinates (see Fig. [I]) is consistent with an overall popu- 
lation of infalling clouds that reflect the sense of Galactic 
rotation, with some peculiar velocities. Although this is not 
a unique interpretation, our finding that the HVC popula- 
tion drops with decreasing z provides additional support to 
this conclusion as this drop is naturally explained if HVCs 



fall through the Milky Way halo onto the disk (see i 3.2 1. A 
better characterization of the neutral and ionized HVC and 
IVC populations as a function of z will be needed to know 
if HVCs actually reach the Galactic disk, but our findings 
further support that they most likely feed the Galactic disk 
with new gas. 

Because they both cover a large fraction of the sky 
and are at distances 5-15 kpc of sun, the HVCs with 
H I column density logA''(H I) <; 18.5 represent a large 
mass reservoir. The total mass of the predominantly ion- 
ized HVCs is M « 1.1 X lO**(d/12kpc)^(/c/O.5)(Z/O.2Z0)-i 
Mq, i.e., about (0.2-2) x 10**(Z/0.2Zq)-^ Mq (see LHll). 
Even though the metallicity can vary from HVC to HVC 
(see above), this implies that the ionized gas mass is at least 
as important - and likely much more important than - the 
H I HVCs with log7V(H I) > 18.5 (excluding the Magellanic 



Strea m), whi ch total about 10 Mq ([W akker 2001] [Putman 
|2006( Wakker et al.|2007| |2008[ [Thorn ct al._2008f . The im- 



portance of the ionized HVCs to the total mass budget is 
also independently demonstrated by the deep Ha emission 



5.3 HVCs in the Local Group and beyond 

Most of the HVCs identifled so far in the Local Group are 
found near galaxies. In the Milky Way direct and indirect 
distance constraints place most of the HVCs and VHVCs 
within about 5-15 kpc (this paper, LHll Wakker|2001 Put 
[man et al.|2003| , with only the Magellanic Stream possibly 
being much farther (Besla et al. 20121. Around M31 and 
M33, the H I HVCs are found within about 50 kpc, and 
many are within 15 kpc ([Thilker et al.||2004| [Westmeier et 



ar] |2005[ |2008[ [Putman et al.||2009[ ). Deep H I observations 
have also revealed a tenuous filament of gas between M31 
and M33 separated by about 150 kpc, which could be the 
result from a tidal interaction between these two galaxies 



(Braun & Thilker 20041. The HVCs observed toward the 



LMC probe outflows from the LMC and are also likely near 
the host galaxy (Lehner et al. 20091. HVCs in the Local 



Group of galaxies, including our Milky Way, are thus all 
within 10-50 kpc. While previously it could be argued that 



purely ionized HVCs could be farther away (e.g., Westmeier 



et al. 20081, our findings show that the ionized and neu- 
tral gas in HVCs are at similar distances, at least in the 
Milky Way halo. Although the sample of galaxies is small, 
there is therefore no evidence in our local galactic neighbor 
that ionized or neutral HVCs are found much beyond 1-2 
galactic radii, except possibly for the remnants of galaxy 
interactions IfJ 

If we move out from our local neighborhood, attempts 
have been made to connect QSO absorbers - in particu- 
larly the optically thick absorbers (a.k.a., the Lyman limit 
systems, LLSs) - w ith the HVCs (e.g., [Charlton et al.|2000 



Richter et al. 2009 1. While some of the LLSs are likely higher 



redshift analogs of HVCs, we emphasize that the distances 
of HVCs place them close to galaxies. HVCs also probe "ab- 
sorbers" with virtually any A''(H I) and therefore likely in- 
clude the population of both weak and strong Mg II ab- 
sorbers. 



^ We also emphasize that while in the Milky Way it might be 
plausible that a population of HVCs may be missed if these have 
low velocities, this is not the case for M31 or M33 where "HVCs" 
with any velocities (included systemic velocities of M31 and M33) 
can be clearly spatially separated from the galaxy. 
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6 SUMMARY 

In this work we have built on the results presented by LHll 
and reported on the covering factors and distances of the 
HVCs, as well as the associations of HVCs observed in H I 
21-cm emission with those found in absorption. We empha- 
size that the UV absorption diagnostics are much more sen- 
sitive to low column density gas than H I emission (or opti- 
cal absorption diagnostics such as Ca II; e.g., Richter et al. 
2011). This allows us to probe clouds over several orders of 
magnitude in H I column density, from A^(H I) < 10^^ cm~^ 
to A'^(H I) > 10'^" cm~'^. We are also able to probe a wide 
range of ionization conditions and densities because we use 
several species in different ionization stages (e.g., O I, C II, 
C IV, Si II, Si III, Si IV). We summarize the results of our 
study Eis follows: 

(i) With a carefully-selected, large sample of AGN that 
cover the sky at |6| > 20°, we find that the covering factor 
of the HVCs with |ulsr| 1? 90 kms~^ of the entire Galactic 
sky is 68% ± 4%. About 74% of the HVC directions have 
iV(H I) < 3 X 10^** cm'2 and 46% have iV(H I) < 8 x 
lO'^'' cm~^ from a comparison of our results with 21-cm H I 
emission surveys. 

(ii) The full sample of AGN and a sample with more uni- 
form and higher S/N data gives the same covering factor. 
HVCs with a total H column density iV(H) < 1Q^\Z/Zq)-'^ 
cm~^ are unlikely to play an important role in the total mass 
of the HVCs. However, HVCs with iV(H I) < lO^'^ cm'^ are 
not rare and have still significant mass in view of the detec- 
tions of metal ion absorption. 

(iii) With the distance, position, and velocity informa- 
tion, we unambiguously associate the denser (H I emission 
selected) and more diffuse (UV absorption selected) regions 
of HVCs. We argue that most of the predominantly neutral 
HVC complexes have ionized envelopes that extend beyond 
the H I contours and that the predominantly ionized HVCs 
contain at least as much mass as the H I component. How- 
ever, there are also large regions of Galactic sky that are 
filled with ionized high-velocity gas with little evidence for 
H I counterparts nearby. One such region is found toward 
the Galactic center and may be generated by Galactic out- 
fiows. 

(iv) The covering factors of HVCs with 90 sC |ulsr| <: 170 
kms~^ determined from the AGN and LHll stellar samples 
are very similar. This confirms that HVCs are within 5-15 
kpc of the sun. The HVCs are therefore fiows of gas in the 
inner Milky Way halo. Our new results also show that the 
covering factor of HVCs drops with decreasing z, which is 
consistent with theoretical predictions of HVCs falling onto 
the Milky Way disk. The HVCs are far enough to have a 
substantial mass (M ex d^), but are also near enough to be 
able to reach the disk. They are therefore the most likely 
source of gas required to support continued star formation 
in the Milky Way. 

(v) While there is no large disparity in the covering fac- 
tor of the HVCs between the two Galactic hemispheres, 
the situation is very different for the VHVCs (HVCs with 
|wlsr| > 170 kms"^). At 6 > 0°, the covering of the VHVC 
is only 13%, but at 6 < 0°, 46% of the Galactic sky is cov- 
ered by VHVCs. Most of the VHVCs are associated with the 
Magellanic Stream at 6 < 0°; a large fraction also appears to 



be related to the leading arm of the Stream (possibly even 
at higher positive latitude than previously thought). Other 
VHVCs with 170 ^ |vlsr| <, 210 kms"^ are afiihated with 
complex A, C, and the anti-center complex. 

(vi) All these elements strongly suggest that in our Local 
Group of galaxies there is no evidence that HVCs are found 
much beyond 50 kpc from their host galaxy, expect for struc- 
tures that are directly related to galaxy interaction (like the 
Magellanic Stream). However, compact VHVCs may have 
been mostly missed along pencil-beam sightlines; these com- 
pact HVCs are the only remaining HVC candidates for dark 
matter mini-halos. 
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APPENDIX A: DISTANCE OF HVC 
COMPLEXES 

We discuss below the connection or absence of connection 
between the HVCs seen in absorption toward the stars from 
the LHll sample (and one additional star HD 195455) and 
in H I emission. 

- NGC5904-ZNG1: No high- velocity H I emission is ob- 
served along this sightline, which is not surprising as |Zech et] 
al. ( 2008 1 derived log A^(H I) = 16.5 using the Lyman series. 
The HVCs at «lsr = —140, —120 kms~^ have similar veloci- 
ties as those observed in complex L seen several degrees away 
(for more information and the spectra, see Zech et al.|2008 |. 
Complex L and this star are in the general direction of the 
Galactic center where lie also the stars NGC6723-III60 (neg- 
ative latitude) and PG1323-086 (positive latitude). Toward 
these stars, high-velocity absorption is observed at about 
—90 kms~^ (see Fig.Mand Table|5|, but they are not close 
to any known H I HVC complex. Similarly, AGN and deep 
H I observations also show additional HVCs at both negative 



and positive velocities toward the Galactic center (Keeney 
et al.|2006 Lockman et al.|2002 1 . It would be interesting to 



map out the entire region of the Galaxy center with deep 
H I and Ha observations at |vlsr| > 90 kms"'^ to deter- 
mine the extent of the ionized gas in this region and its 
possible connection with complex L. The HVCs observed 
toward NGC5904-ZNG1 have a very well determined metal- 
licity based on O I and H I absorption, which is super-solar 
([0/H] = +0.22 ±0.10). (Another test for the connection or 
not of these clouds with complex L would be to measure the 
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metallicity of complex L). In similar directions, O VI absorp- 
tion is also detected at these positive and negative velocities 



closer to the plane (e.g., toward HD168941, see Bowen et al, 
2008), so these HVCs could be affiliated with flows originat- 



ing in the disk. The super-solar metallicity, large column of 
predominantly ionized gas seen in the direction of the inner 
Galaxy, presence of both positive and negative velocity may 
suggest that they participate to a galactic fountain fed by 
outflows occurring in the central regions of the Milky Way. 
- M13-Barnard29: An HVC is observed in absor ption 
at -121 km s"^ (see the spectra in [Welsh et al.||201l| ). The 
sightline passes near to (but outside of) the H I contours of 
complex C and no LAB H I emission is observed at —121 
kms~^. However, the absorption velocity is similar to veloc- 
ities observed in the nearby H I emission from complex C. 
The distance of complex C was derived to be at d = 10 ±2.5 
kpc (Thom et al. 2008, and see also Wakker et al. 2007). 
While the distance of M13-Barnard29, 7.1 kpc, is slightly 
outside this distance bracket, it is not inconsistent, espe- 
cially taken into account a 10% error on the distance of M13- 
Barnard29, possible variation in the distances of complex C 
with position, and similarity in the velocities. Therefore this 
sightline probably traces the extended ionized envelope of 
complex C. We therefore disagree with the conclusions of 



Welsh et al. ( 2011 1 who argued that the distance of this star 



was too small for the gas to be associated with complex C. 
- HD195455: This star lies in the general direction of 
the IVC complex gp, which spans LSR velocities between 
-1-55 and -1-90 kms"\ The HVC observed at -f95 kms"^ 
toward this direction could be a higher velocity extended 
ionized component of this complex. The IVC Complex gp 
component is also observed at -|-70 kms~^ along this sight- 



line (see Fig. 111. The distance of HD195455 is consistent 
with the bracketed distance of complex gp between 1.8 and 
3.8 kpc ( |Wakker et aL]|2008[ ). 

- HS1914+7139: Two HVCs were detected along this 
sightline at vlsr = -118,-175 kms~^, which are associ- 
ated with the Outer Arm and/or complex C (see for more 
information and the spectra in |Lehner fc Howkl|2010| and 
also discussion in Tripp & Song 2012). LAB H I emission 
was detected only in the —112 kms~^ component. 

- PG0122+214: This hne of sight has also two HVCs 
detected in absorption at —91 and —160 kms^^ (see Fig.l9| 
and passes near the anti-center HVC complexes, more specif- 
ically near the Cloud WW507 and the Cohen stream, where 
similar velocities are observed. Again no high-velocity H I 
emission is observed along this sightline. The distance of 
PG0122+214 (9.6 kpc) is cons istent with the brack et dis- 
tance of 5-11.7 kpc derived by Wakker et al. (20081. How- 



ever, the higher velocity component was not observed in the 



stellar spectra of Wakker et al. (20081, so this star places a 



new upper limit to the distance of the higher velocity gas 
observed in this complex. 

- PG0832+675: This line of sight diflters with the oth- 
ers in the sense that H I emission is observed with the peak 
emission centered at —145 kms~^, which is different from 



the observed absorption centered at —123 kms . In Fig. 10 
we show the LAB H I 21-cm emission and the normalized 
spectra of O I, C II, Al III, Si II, S II, and Fe II. The central 
velocities of the low and intermediate velocity components 
seen in S II absorption and H I emission are the same, and 
therefore the difference between the absorption and emis- 



sion at high velocities is real. Although the UV spectrum 
is complicated by many stellar features, the component at 
— 123 kms~^ is very likely interstellar because the veloc- 
ity is the same for different species (C II, O I, N I, Si II, 
Fe II, although we note that O I and N I are likely contam- 
inated at —123 kms~^) and different from the star (—70 
kms~^), and because the column densities derived from dif- 
ferent transitions of the same species give the same values 
(logiV(Fe II) = 12.65 ± 0.02 based on the A2344, 2382 tran- 
sitions and logiV(Al III) = 12.55 ±0.05 based on the A1854, 
1862 transition, and logiV(Si II) = 13.00±0.08 based on the 
A1193, 1260, 1304 transitions). 

PG0832±675 lies in the direction of complex A within 
the H I contour, so it is not surprising to see H I emission in 
this direction. The star is likely at the front edge of complex 
A, which would explain why the O I is so weak at —145 
kms~^ (the metallicity of complex A -B.P. Wakker, private 
communication 2011 - would imply a far stronger absorp- 
tion). Wakker (2001 1 summarized a distance bracket for this 
complex of 4.0-9.9 kpc (see Wakker et al.|[T996 van Woer- 



den et al.|[l999[ ). Given that we see little O I, C II, Si II, 
and Fe II absorption at -145 kms"^ toward PG0832±675 
compared to the H I column density at the same velocity, 
we conclude that most of the H I lies behind the star. This 
implies a distance bracket for the predominantly neutral gas 
of complex A of 8.1-9.9 kpc. The —123 kms~^ component 
may be the ionized envelope of complex A or an unrelated 
foreground cloud. The upper limit for its distance, in either 



case, is thus 8.1 kpc. While Barger et al. ( 2012 1 observed Ha 
emission in complex A at LSR velocities ranging from —220 
to —110 kms~^, there is no strong evidence of Ha emission 
at —123 kms"'^ over a 1° beam toward this star (K. Barger 
2011, private communication). This does not rule out that 
the —123 kms~^ component is part of complex A. 

- PG1002+506: The H I HVC WW503 is about 2.5° 
away from this sight line with ulsr = —105 kms~^ (B.P. 
Wakker, private communication 2011), similar to the veloc- 
ity observed in absorption. Our observations provide a new 
upper limit, < 13.9 kpc, on the distance to the ionized com- 
ponent of this cloud. This strongly suggests that most of the 
HVCs in the anti-center region at 6 < 0° are within about 10 



kpc (see Wakker et al. 2008 1. On the other hand, no H I emis- 



sion is observed near this direction at ^lsr = ±101 kms^ , 
indicating the presence of a positive velocity HVC that is 
not associated with a H I complex at a distance d < 13.9 
kpc. 

- PG0855+294: Several small high- velocity H I clouds 
are observed in this general direction. This sightline passes 
0.5° away from HVC WW113 seen in H I emission at 
«LSR ~ —90 kms~^(B.P. Wakker, private communication 
2011), setting an upper limit of 6.5 kpc for the extended 
ionized component of these clouds. 

- PG09 14+001: This sighthne is near complex WA 
(WLSR = ±140-±170 kms~^) and WB (ulsr = ±80-±110 
kms~^), and absorption features at ±110 and ±170 kms~^ 
are observed in the spectrum of this star (for more informa- 
tion and the spectra, see LHll). This places an upper limit 
of 16 kpc for the extended ionized envelope of these HVCs. 
For the WB complex, the HVCs could be closer as [Thorn et| 
al. (20061 derived a bracket distance of 7.5-11. Ikpc. 



- EC10500-1358: This sightline is also in the general 
direction of complex WA, but the velocity seen in absorption 
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(wlsr = +97 kms ^, see the spectra in Fig.l9| is smaller. It 
is, however, about 0.5° from HVC WW95 at «lsr = +100 
kms~^, setting a distance limit of 5.2 kpc for this cloud. 

- NGC104- UIT14: This sightline shows no HVC absorp- 
tion, although it has strong H I 21-cm emission observed at 
+ 112 kms"^ , with a H I column density logiV(H I) = 20.39. 
This sightline passes in the southern part of the WE com- 
plex ( |Morras et~aL]|2000| | Wakker|[200T| ) , which is a patchy 
ensemble of HVCs at similar velocities rather than a large 
continuous structure like complex C. The limit from O I 
gives \ogN{0 I) < 13.8 at 3(7. It is unlikely that the metal- 
licity in this complex is < 0.5 x 10"'^ ^0, implying that the 
star sets a lower limit on the distance of WE, and hence us- 
ing the upper limit from |Sembach, Savage, fc Massa| ( |1991[ ), 
the WE complex is at 4.5 < d < 12.8 kpc. 
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